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ABSTRACT 
The  theory  and  equations  for  forecasting  the  wind  drift  of  sea  ice, 
as  previously  developed  by  several  researchers,  are  being  used  in  the  fore- 
casting of  that  component  of  sea-ice  drift.   The  existing  equations  are 
modified  and  programmed  for  use  with  a  computer  to  produce  a  forecast 
of  wind  drift,  thus  permitting  rapidly  made  forecasts  over  large  areas ■ 
Applicability  and  limitations,  ways  of  handling  other  factors,  i.e..    per- 
manent currents  and  melting  and  freezing,  and  methods  for  general  im- 
provement of  the  programs  are  discussed.   The  programs  developed  are  use- 
ful as  research  tools  and  as  applied  to  forecasting . 
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lo   Introduction. 

The  problem  of  forecasting  the  wind  drift  of  sea  ice  has  been  im- 
portant to  man  since  he  first  attempted  to  navigate  the  frozen  sea.   Use 
of  the  polar  sea  for  supply  and  commercial  routes  and  as  an  operating 
area  for  nuclear  submarines  makes  desirable  faster  methods  of  ice  fore- 
casting for  large  areas.   Zubov,  Shuleikin  and  Fukutomi  [1,  2,  3],  among 
others ,  have  investigated  the  problem  of  wind  drift  of  sea  ice  and  have 
proposed  various  formulae  for  prediction  of  its  movement.   The  U.  S.  Naval 
Oceanographic  Office  (formerly  U.  S.  N.  Hydrographic  Office)  has  published 
ice-forecasting  techniques,  including  those  for  forecasting  wind  drift, 
which  are  used  by  Navy  meteorologists  and  others  [4]. 

The  Oceanographic  Office  method  employs  graphs  and  tables  to  obtain 
a  forecast  by  hand  calculation.   Forecasting  for  any  extensive  area  is  a 
slow  process  and  much  subjectivity  is  involved.   Data  available  for  hand 
forecasts  are  rather  sparse,  making  it  difficult  to  do  a  complete  forecast 
even  if  time  were  not  a  factor. 

A  machine  method  of  predicting  the  wind  drift  of  sea  ice  would  represent 
a  significant  improvement  over  present  techniques;  (1)  it  would  rapidly 
provide  a  complete  wind-drift  forecast  for  an  area;  (2)  continuity  of  fore- 
casts would  result  from  retention  of  computed  values  which  may  be  modified, 
in  part,  by  new  observations  and  carried  forward  into  new  forecasts;  and 
(3)  forecast  wind  information,  to  serve  as  input  data,  is  available  via 
numerical  methods. 

This  paper  will  present  one  approach  to  the  numerical  forecasting  of 
wind  drift  of  sea  ice,  together  with  an  examination  of  its  applicability 
and  limitations. 


2 „    Background. 

Many  expeditions  to  the  Arctic  have  contributed  to  available  infor- 
mation on  the  wind  drift  of  ice,  geology,  meteorology  and  oceanography 
of  this  area;  these  data  have  been  used  by  later  researchers.   In  September 
of  1893  Nansen's  ship  FRAM,  built  to  endure  the  crushing  pressures  of  ice 
floes,  entered  the  pack-ice  in  latitude  78-50  N  off  the  New  Siberian  Islands 
for  the  express  purpose  of  gathering  data  on  the  motion  of  ice  and  the  currents 
In  August  of  1896  after  drifting  for  three  years  with  the  ice  and  reaching 
85-55  N,  farthest  north  for  any  surface  ship,  FRAM  emerged  from  the  ice 
into  the  Norwegian  Sea  east  of  Greenland  completing  its  gathering  of 
scientific  data  in  the  Arctic.   From  1937  to  1940  the  Russian  icebreaker 
SEDOV  drifted  over  a  route  similar  to  that  of  FRAM„   In  the  intervening 
years  several  vessels  had  completed  shorter  drifts. 

During  1937  a  four -man  Russian  team  under  Ivan  Papanin  drifted  for 
nine  months  on  an  ice  floe,  the  Papanin  Floe;  thus,  the  Russians  pioneered 
in  placing  scientific  stations  on  drifting  ice  islands  [Dyson,  5].   Since 
1950  a  number  of  these  drifting  stations  have  been  occupied  by  the  United 
States  and  by  Russia. 

Zubov  [1]  utilized  data  obtained  by  the  SEDOV  expedition  to  obtain 
empirical  rules  for  the  wind  drift  of  close  ice   (ice  which  covers  most 
of  the  surface).   His  results  compared  remarkably  well  with  those  which 
Nansen  derived  from  data  obtained  during  the  drift  of  FRAM,   Zubov 
concluded  that  over  a  deep  sea,  away  from  the  influence  of  land,  the  wind- 
caused  ice  drift  is,  on  the  average,  at  an  angle  of  29   to  the  right  of  the 
wind  and  at  a  speed  which  is  1/50  of  the  speed  of  the  wind.   In  his  treatment 
of  the  drift  of  open  ice  (ice  covering  between  two  and  six  tenths  of  the 
surface)  Zubov  shows  that  the  speed  of  drift  is  a  function  of  concentration 

and  increases  with  both  decreasing  concentration  and  increasing  roughness. 
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Shuleikin  [6],  using  dynamic  relations  and  the  drift  of  the  Papanin 
Floe  compared  to  winds  derived  a  theory  of  ice  drift  in  the  open  ocean. 
In  a  later  writing  [2]  he  examines  his  theory  further  by  using  data  from 
the  Papanin  and  NP-1  floes  and  obtains  constants  for  his  equations. 
Shuleikin's  comparatively  simple  equations  consider  only  ice  fields  of  10/10 
concentration. 

Since  Zubov  does  not  give  equations  for  the  speed  of  ice  of  varying 
concentration,  we  refer  to  Fukutomi  [3]  for  a  treatment  of  wind  drift  of 
sea  ice  under  these  conditions.   Fukutomi' s  formulae  are  considerably 
more  complicated  than  those  of  Shuleikin,  introducing  stresses  on  the 
vertical  surfaces  of  floes  by  wind  and  current  and  stress  of  wind  on  the 
exposed  water;  Shuleikin  considered  only  stresses  on  the  horizontal  ice 
surfaces.   This  treatment  by  Fukutomi  makes  necessary  auxiliary  equations 
and  more  data  on  floe  dimensions .   Further,  Fukutomi's  constants  were 
determined  from  data  obtained  in  the  Sea  of  Okhotsk  and  thus  differ  from 
Shuleikin's  constants,  based  on  Arctic  Ocean  data.   However,  similarities 
exist  in  the  theories  and  results.   Fukutomi  concludes  that  movement 
predicted  by  his  equations  agrees  well  with  Nansen's  observations.   His 
results  support  Zubov  in  that  drift  speed  increases  with  decreasing  con- 
centration.  Comparison  of  Fukutomi's  results  for  10/10  concentration  with 
those  of  Shuleikin  shows  rather  good  agreement,  particularly  for  the  angle 
of  deviation  of  ice  drift  from  wind  direction   Deviation  angle  is  also 
shown  to  depend  on  concentration;  this  is  a  reasonable  result,  since  an 
increased  wind-induced  current  in  exposed  water  areas  directly  affects 
ice  movement. 

Review  of  the  above  papers  and  others  indicates  that  the  several 

theories  and  methods  are  mutually  consistent;  their  different  results  are 

related  largely  to  local  differences  in  ice  characteristics  (e.g.  roughness) 
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which  are  described  by  certain  parameters.   Therefore,  the  method  developed 
here  uses  the  equations  of  Shuleikin  modified  for  concentration  by  reference 
to  Fukutomio 


3.   Models  for  the  wind  drift  and  changes  in  concentration  of  sea  ice. 

Shuleikin  [2]  has  provided  the  basic  equations  required  to  obtain 
close  ice  movement.   These  are: 

tan  if   =  tan  8  +  B/u  (1)  , 


and   u  =  V 


(§  cos  a  +  cos  B) 


k  6 
a  a 

k  6 
w  w 


cos  y 
cos  8 


(2) 


where 


if   -   angle  between  wind  direction  +180  and  the  direction  of 
ice  drift 

u  =  ice  drift  speed 

V  =  wind  speed 

k  =  coefficient  of  friction  between  air  and  water 
a 

k  ■  coefficient  of  friction  between  ice  and  water 
w 

6  =  density  of  air 

a  J 

6   =  density  of  water 
w 

co  =  angular  velocity   of   the   earth 
0  =   latitude 


co  =  co  sin 


A  =  a  dimensionless   coefficient   =   1.27   x   10 


N  =     rr  k     6 
a     a 

nr-   .2 


V  2  to  A     6 


w 


5-    nk 


w 


N  Id  JT 
a  =  actn  (§N/T+   1) 
8  ■  rr/4   -  a 

m  =  mass    of   ice   per  unit   surface   area 

_  2 

B  =  2m  cd   (5   cos   a  +  cos   g) 

k      6      cos    ft 

WW 


Figure  1  depicts  the  relationships  among  some  of  the  above  quantities. 
In  deriving  the  above  equations,  Sheleikin  has  assumed  an  ice  concentration 
of  10/l0o   He  further  assumed  the  ice  drift  is  in  equilibrium  with  a  constant 
wind,  i.e.  an  unaccelerated  motion  prevails,  and  all  stress  is  applied  to 
the  horizontal  surfaces  of  the  ice. 

Fukutomi's  equations  require  a  knowledge  of:   (1)  total  area  of  ice, 
(2)  area  of  exposed  sea,  (3)  plane  surface  area  for  each  floe,  (4)  vertical 
area  exposed  to  wind  pressure  for  each  floe,  and  (5)  vertical  area  exposed 
to  current  for  each  floe*   Moreover,  one  must  be  able  to  obtain  or  estimate 
the  energy  loss  due  to  eddy  currents.   Most  of  these  factors  are  not  known 
by  observation  and  are  difficult  to  estimate  with  any  degree  of  accuracy; 
therefore,  it  is  not  practicable  to  use  Fukutomi's  equations  idrectly„ 
The  equations  are  too  long  to  reproduce  here,  and  the  reader  is  referred 
to  Fukutomi  [3]  for  further  study. 

It  is  possible,  however  to  approximate  one  set  of  Fukutomi's  curves 
for  wind  factor,  u/V,  versus  wind  speed,  V,  with  Shuleikin's  equations  by 


1  +  (  "^p  )  U  +  3V'  x  >01] 


the  introduction  of  a  correction  factor,  E  - 

where  c  is  ice  concentration  and  V  is  wind  speed  in  knots,  (c  <  10/10). 

The  approximation  is  shown  in  figure  2  where  the  curves  according  to  Fukutomi 

represent  complete  energy  loss  due  to  eddy  currents;  this  corresponds  well 

to  that  condition  when  area  of  ice  is  large  compared  to  ice  thickness.   The 

resulting  equations  are : 

tan  if   =   tan  0  +  B/uE  (3), 


T,r/s       .     «\   a  a  cos  * 

and   u  =  V[(§  cos  a  +  cos  fl)     : 7—  f- 

1 N^  k  0  cos  6  J 

w  w 


(4). 


The  modification  was  obtained  in  a  purely  arbitrary  manner  in  order  to  best 
fit  Fukutomi's  curves  over  the  range  of  winds  normally  observed. 
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Now,  a  method  is  desired  to  obtain  the  rate  of  change  of  ice  con- 
centration at  a  point  and  the  resulting  concentration  at  that  point  for 
some  future  time.   Simpson  [7]  provides  equations,  which  are  readily 
adaptable  to  numerical  methods,  for  computing  these  quantities.   The 
equation  of  continuity  of  ice  concentration  given  by  Simpson  is: 

dc  _    /   i  .  L  \      cte      dc  /CN 

¥  ■c(^  +  ^)"ui^"v15  (5)> 

and  the  equation  for  concentration  at  a  given  time  is ; 

C(v  Vj+i>  =  c<*i.yi*j>  +  cl?(xi'yi'¥j)1,[tj+rtj1  (6); 

where : 

dc 

t-  =  local  change  in  concentration  with  time 

c  =  concentration 

u.=  x-component  of  ice-drift  velocity 

v.-  y-component  of  ice-drift  velocity 

x  =  x-coordinate  of  location 

1 

y.=  y-coordmate  of  location 

t.=  initial  time 
J 

t .  ,  .,  =  time  at  which  new  concentration  is  computed 

J+l 

be 
t.  ■  time  over  which  — ■  is  computed,  the  time  step. 

J  at 

Equations  (3),  (4),  (5)  and  (6)  permit  one  to  determine  wind  drift 
of  sea  ice  and  concentration  changes  resulting  from  that  drift,  thus 
completing  our  models. 


4.      Development   of   the   computational   programs. 

One   of   the    two   programs,    ICEF0R1 ,   might   be   considered  as   a   step   to- 
ward  the   final   program  developed  here,    using   the  models   developed    in 
section   3.      ICEF0R1 ,    reproduced    in  Appendix  A,   was  written  using   equations 
(1)   and    (2)         (Shuleikin's   equations)  ;    it   produces  wind-drift   vectors 

for   ice    fields  with  k     and  k     known  and   uses    inputs    of    latitude   and  wind 

w      a  r 

velocity.   While  it  is  not  the  final  form  of  the  wind-drift  forecast,  it 
does  have  some  utility. 

The  next  step  was  to  program  equations  (3)  and  (4),  those  modified 
for  concentration,  to  obtain  outputs  of  ice  drift  similar  to  the  output 
from  ICEF0R1 .   Basically  a  program  segment  of  ICEF0R1  was  modified  to 
accept  all  ice  concentrations.   Changes  in  concentration  were  next  con- 
sidered.  Equation  (5)  is  put  into  finite  difference  form: 

Ac         ruru3  +  v2-v4,       .  crc3,         rVC4  . 

At"  =   "Co[  ~W  +  ~^D~]    "   Uo[  ~ 2D~]    "   Vo    [~2b~]  (7)> 

where : 

AC  .  r  r- 

—  =  time  rate  of  change  of  concentration 
At  e 

c   =  concentration 

u  =  x- component  of  drift  velocity 

v   =  y-component  of  drift  velocity 

D  =  distance  between  data  points  and  subscripts  indicate  data 
points  as  follows : 


1.     .     .3 


Solution  of  equation  (7)  requires  values  of  ice-drift  velocity  (u  ,v  ) 

and  initial  concentration,  c  ,  at  each  data  point  for  each  step.   A  new 

10 


concentration,  c  +  Ac,  can  then  be  determined  for  each  data  point.   One 
o 

Ac 
notes,  however,  that  in  order  to  compute  -rr  at  any  one  point9data  from 

four  other  points  are  required.   To  satisfy  this  requirement,  and  to  obtain 

an  area  forecast,  a  rectangular  grid  to  approximate  an  area  of  the  earth's 

surface  was  devised.   The  grid  is  a  546-point  array  of  dimensions  21x26. 

It  has  an  adjustable  scale  which  permits  the  distance  between  grid  points 

to  be  changed  by  use  of  an  appropriate  scale  factor.   Obviously  if  this 

plane  is  made  too  large  it  no  longer  approximates  an  area  of  the  earth's 

spherical  surface.   However,  for  an  area  the  size  of  Baffin  Bay,  requiring 

a  scale  factor  of  2.0  for  representing  30  nautical  miles  between  grid 

points,  the  approximation  to  the  earth's  surface  is  good  and  no  problems 

arise.   The  grid  permits  input  of  data  over  an  area  and  further  permits 

computation  of  convergence  and  advection  within  that  area,  as  discussed 

in  detail  later.   The  program  is  set  up  to  determine  ice-drift  vectors, 

changes  in  concentration  and  new  concentration;,  point  -for-point  throughout 

the  grid;  it  is  named  ICDGRID. 

Several  authors,  notably  Browne  and  Crary  [9]  and  Fukutomi  [3], 

have  obtained  values  for  the  coefficients  of  friction  (k  and  k   in 

a      w 

equations  (2)  and  (4))  which  differ  from  those  used  by  Shuleikin,   Also, 
Wittmann  and  MacDowell  [4]  apparently  use  different  coefficients  in  modi- 
fying Shuleikm's  equations  for  Baffin  Bay.   It  is  apparent  that  these 
coefficients  do  vary,  as  we  might  expect,  for  there  are  variations  in  the 
roughness  of  ice  with  age  (newly- formed,  thin  ice  is  smoother  than  old 
ridged  ice)  and  with  location.   Therefore,  the  programs  developed  here 
permit  the  user  to  enter  values  of  these  coefficients  which  best  fit 
the  area  of  forecast. 

Special  problems  arise  when  land  boundaries  are  approached ,"  therefore 
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the  program  is  written  to  recognize  points  of  concentration  entered  as 
2.0  as  land  points  and  to  ignore  these  points  in  computing  ice  drift. 
Certain  of  the  land  points  are  required s  however,  in  the  finite  differencing 
process  used  to  obtain  cWc3t„   In  order  to  treat  concentrations  near  the 
boundaries  the  program  is  written  to  see  the  concentration  of  land  points 
as  10/10  when  the  drift  is  onshore  and  as  0/10  when  the  drift  has  an 
offshore  component.   In  this  manner  changes  in  concentration  along  land 
boundaries  are  made  realistic . 

Program  ICEGRID  is  outlined  in  Appendix  B.   It  utilizes  inputs  of 
ice  concentration,  ice  thickness  and  mean  daily  surface-wind  velocity  to 
obtain  a  five-day  wind-drift  ice  forecast  in  one-day  steps.   The  outputs 
are  grids  of  ice  concentration,  direction  of  movement,  and  24-hour  distance 
of  movement  for  each  of  the  five  days  of  the  forecast. 
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5.  Applicability 

The  programs  developed  here  are  applicable  to  any  geographical  area. 

If  properly  indexed  in  latitude,  ICEGRID  will  compute  over  the  geographical 

area  desired,  and  ICEF0R1  computes  for  geographical  points  as  directed. 

Accurate  forecasting  in  each  case  requires  that  the  coefficients  of 

friction,  k  and  k  ,  which  tend  to  vary  with  location,  be  known.   Then, 
'aw 

given  primary  inputs  of  ice  concentration,  ice  thickness  and  wind  velocity, 
proper  outputs  may  be  produced. 

There  are  two  basic  uses  of  the  programs :   forecasting  the  wind-drift 
component  of  sea-ice  drift  and  research.   Both  ICEGRID  and  ICEF0R1  provide 
forecasts  of  the  wind  drift  of  sea  ice,  with  ICEGRID  providing  the  most 
complete  forecast.   ICDGRID  provides  an  area  forecast  and  gives,  in  addition 
to  direction  and  distance  of  movement,  a  concentration  field:  ICEF0R1  produces 
only  the  ice-drift  steering  vectors  at  selected  geographical  points. 

With  regard  to  research  ICEF0R1  may  be  used  to  obtain  paired  values 

of  k  and  k   for  an  area  in  which  sufficient  observational  data  of  actual 
a      w 

drift  are  available  for  comparison  with  computations.   The  rapid  computa- 
tional ability  available  in  machine  methods  makes  this  an  extremely  advantageous 
tool.   Additionally,  one  may  easily  extend  the  possibilities  for  research 
with  both  programs  to  include  investigation  of  permanent  currents  in 
the  frozen  sea,  to  study  effects  of  specific  features  (i.e.  pressure  ridges, 
floe  size)  on  the  wind  drift,  and  to  relate  drift  to  geographical  features. 
The  fact  that  ICEGRID  is  applicable  to  all  concentrations  permits  its  use 
in  more  general  research  than  ICEF0R1,  which  is  restricted  to  concentra- 
tions of  10/10. 
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6.   Limitations. 

This  section  discusses  the  limitations  of  ICEGRID  as  a  forecasting 
method „   Since  only  a  limited  portion  of  the  total  forecasting  problem 
is  treated  in  the  theory  used  here,  only  that  part  of  the  drift  directly 
due  to  the  wind  is  found;  components  due  to  permanent  currents  and  concen- 
tration changes  due  to  melting  and  freezing  are  ignored „   Some  of  the 
limitations  discussed  are,  therefore,  not  related  to  the  stated  problem, 
that  of  forecasting  the  wind  drift  of  sea  ice. 

Permanent  currents. 

Probably  the  most  important  limitation  of  the  program  is  the  failure 
to  consider  the  effects  of  permanent  currents.   Wherever  such  currents 
exist  they  must  affect  the  total  drift  of  the  ice „   Therefore,  the  program 
should  be  expanded  to  include  these  effects. 

One  may,  of  course,  subjectively  superimpose  these  effects  upon  the 
computer  output,  drawing  upon  a  knowledge  of  the  permanent  currents  to 
infer  the  total  ice  movement .   A  more  direct  method  of  considering  these 
currents  would  be  to  input  permanent  current  vectors  at  each  grid  point, 
adding  them  vectorially  to  the  ice  drift  due  to  wind.   The  total  ice 
drift  and  the  resulting  changes  in  concentration  would  then  be  given  by 
the  computer  output. 

Curiously,  permanent  current  data  are  not  always  readily  available. 
Most  current  observations  are  made  during  periods  when  little  or  no  ice 
is  present,  and  these  would  likely  be  poor  approximations  for  times  of 
the  year  when  considerable  ice  cover  is  present.   Still,  in  many  areas, 
say  in  the  East  Greenland  Current  or  in  the  Labrador  Current  area,  a 
calculation  based  on  any  approximation  of  the  permanent  currents  would 
surely  be  of  more  use  than  one  considering  wind-drift  alone. 
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Melting  and  freezing. 

A  problem  perhaps  equal  to  that  of  permanent  currents  is  that  of 
assessing  the  influence  of  melting  or  freezing  during  the  forecast  period. 
Assuming  no  effect  may  be  reasonable  early  in  the  melting  season;  however, 
it  is  most  unrealistic  at  the  height  of  the  melting  season  and  early  in 
the  freezing  season.   At  these  times,  ice  boundaries  may  actually  move 
opposite  to  the  drift  if  these  effects  are  large  enough.   This  time  of 
year  from  greatest  melting  to  early  freezing  is  also  the  most  important 
to  the  forecaster?  as  this  is  the  time  when  ships  ply  the  ice  covered 
sea  and  the  need  for  an  accurate  forecast  is  greatest.   Although  the 
problems  are  complex,  models  which  account  for  ice  melting  and  freezing 
exist  and  the  necessary  meteorological  and  oceanographic  data  may  be 
extrapolated  or  approximated  from  climatology.  A  machine  method  of  ac- 
cumulating degree  days  of  frost  or  warming  degree  days  might  then  lead  to 
machine-produced  estimates  of  freezing  or  melting  to  be  incorporated  into 
the  ice  forecast;  thus  a  more  complete  forecast  would  result. 

Hummocking  and  ridging. 

The  effects  of  additional  hummocking  and  ridging  are  neglected  once 

k  and  k  are  chose;  actually,  the  drift  characteristics  of  the  ice  are 
a      w 

altered  through  changes  in  the  roughness  parameters.   The  immediate  re- 
sult of  failing  to  consider  the  above  factors  is  a  forecast  of  slower 
than  actual  drift  in  areas  characterized  by  increased  roughness;  for, 
as  previously  stated,  Zubov  [1]  concludes  that  increased  hummocking  increases 
drift  speedo   In  Fukutomi's  development  [3]  hummocking  and  ridging  are  con- 
sidered in  that  they  increase  vertical  areas  exposed  to  both  wind  pressure 
and  currents   Clearly,  these  factors  should  be  given  some  consideration; 
but  with  presently  available  data  any  modification  would  be  largely  sub- 
jective, as  aerial  observations  indicate  only  the  presence  and  not  the 
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dimensions  of  such  features;  further,  one  may  not  know  the  extent  of  these 

features  present  in  the  data  from  which  the  determination  of  k  and  k 
r  aw 

was  made . 

It  is  also  highly  desirable  to  be  able  to  predict  the  formation  of 
ridges  and  hummocks,  as  ice  dominated  by  these  forms  is  difficult,  and  per- 
haps impossible,  to  navigate.   Large  hummocks  and  ridges  also  pose  a  hazard 
to  submarine  navigation  by  restricting  water  depth.   They  may  lend  to  in- 
creased formation  of  fast  ice,  if  they  ground  in  shallow  water,  and  thus 
alter  the  effective  land  boundary «   This  problem  requires  further  study 
directed  toward  a  dynamical  solution,  but  will  not  be  considered  here „ 

Transport  from  outside  the  area. 

No  consideration  is  made  of  ice  advected  into  the  forecast  area  from 
beyond  its  boundaries,  since  no  means  is  provided  for  entering  data  from 
outside  the  area  into  the  machine .   This  inability  to  compute  new  con- 
centrations at  the  grid  boundary  points  leads  to  corresponding  errors  in 
concentration  near  the  boundaries  .   The  problem  may  be  further  complicated 
by  iterated  use  of  the  output  as  continuity  data  for  new  forecasts.   If 
boundary  points  are  not  modified  by  new  observational  data  anomalous  fore- 
casts near  the  boundary  can  result,  since  the  original  input  concentrations 
will  have  been  retained  at  the  boundary  points.   Accordingly,  the  numerical 
forecaster  must  be  alert  to  change  such  points  from  climatology  when  such 
changes  are  required  in  the  interest  of  continuity  and  common  sense.   Of 
course,  one  might  use  a  grid  scale  to  forecast  over  an  area  larger  than 
that  actually  required;  thus,  the  boundary  points  would  have  less  effect 
on  the  desired  forecast  area. 

Nearshore  drift. 

In  the  nearshore  region  of  a  water  body,  including  that  portion  near 
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the  edge  of  the  fast  ice,  the  drift  of  sea  ice  may  be  quite  unlike  that 
in  the  open  ocean.   Zubov  [1]  discusses  the  nearshore  drift  with  respect 
to  the  effect  of  winds  in  piling  up  or  driving  away  the  ice  from  the  shores. 
Shuleikin  compares  ice  drift  along  East  Greenland  with  that  in  the  open 
Arctic  and  shows  the  differences  which  exist.   From  these  and  other  dis- 
cussions, it  is  apparent  that  the  presence  of  a  solid  boundary  causes 
local  changes  in  concentration  and  alters  the  direction  and  speed  of  ice 
drift. 

The  effect  of  a  solid  boundary  on  concentration  is  to  act  as  a  nucleus 
for  increasing  concentration  when  there  are  onshore  winds;  and,  contrarily, 
to  act  as  a  zero  concentration  point  in  the  case  of  offshore  winds  (there 
is  no  source  of  sea  ice  on  the  land  to  contribute  to  offshore  advection). 
This  portion  of  the  problem  is  considered  in  ICEGRID  by  forcing  the  computer 
to  see  concentrations  at  adjacent  land  points  as  0/10  and  10/10  for  offshore 
and  onshore  winds,  respectively,  in  computing  changes  in  concentration  at 
nearshore  grid  points. 

The  influence  of  the  boundaries  on  the  drift  vector  itself  is  not  con- 
sidered,  To  make  such  consideration  would -require  that  the  program  have 
available  not  only  the  location  of  land  points  but  the  orientation  of  the 
boundary  as  well.   To  insert  such  information  into  the  program  would  have 
required  tying  it  to  a  specific  geographical  location  which  would  have  made 
it  less  general;  therefore,  no  attempt  was  made  to  handle  this  effect. 

Thickness . 

The  input  of  ice  thickness  is  as  much  a  problem  as  considering  hummocks 
quantitatively.   Few  observations  of  thickness  are  available,  from  an  oc- 
casional ship  report  during  the  navigable  season  or  from  a  land  station; 
these  stations  may  report  regularly  on  fast  ice  in  the  area.   The  result 
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is  that  thickness  must  be  input  from  climatology  or  long-range  ice  outlooks 
formed  from  climatology  and  heat  budget  considerations.   These  provide 
reasonably  good  data  over  large  areas5  but  local  anomalies  may  be  entirely 
smoothed  out,  especially  by  the  use  of  climatology .   Fortunately  my  compu- 
tations using  Shuleikin's  equations  show  that  drift  distance  decreases 
on  the  order  of  0.5  nautical  miles  per  day  as  ice  thickness  increases  from  5 
to  8  feet s  for  all  wind  speeds.   The  angle  of  deviation  from  the  wind  increas- 
es with  increasing  thickness;  however,  for  a  thickness  change  from  5  to  8  feet 
the  maximum  difference  in  the  angle  is  about  15  at  low  wind  speeds  and  less 
with  greater  wind  speeds,  being  on  the  order  of  5   for  20-knot  winds.   Thus, 
the  climatology  or  ice  outlook  thicknesses  will  give  generally  acceptable 
results  when  used  in  computations. 

Transient  states  in  the  ice  field. 

In  the  equations  for  sea  ice  drift  it  is  assumed  that  the  ice  moves 
everywhere  in  equilibrium  with  the  mean  wind  over  the  period  of  a  forecast 
step,  in  this  instance  one  day.   It  should  be  apparent  that,  due  to  inertia 
of  the  ice  floes,  a  transient  stage  occurs  whenever  a  change  in  wind  velocity 

occurs.   Zubov,  however,  states  that,  "all  observers  are  struck  by  the  speed 

1 
with  which  ice  reacts  to  any  change  in  the  wind,"   indicating  transient 

states  in  the  ice  movement  to  be  rather  short.   They  are  ignored,  therefore, 

since  the  program  steps  are  of  one-day  duration.   It  might  be  necessary  to 

consider  these  states  for  shorter  program  steps. 

Observational  data. 

Any  forecasting  program  may  be  severely  limited  by  the  availability, 
accuracy  and  form  of  data  to  be  used  as  input.   Currently  available  obser- 
vations of  sea  ice  limit  the  forecast  in  all  of  these  respects.   The  majority 


Zubov,  N.  N. ,  Arctic  Ice.   Moscow,  1953.  p.  355, 
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of  data  are  from  aerial  reconaissance  flights  made  at  somewhat  irregular 
intervals  over  the  limits  of  the  ice;  thus  the  amount  of  data  available 
is  adequate  over  only  a  small  area. 

In  addition  to  data  scarcity,  data  reliability  is  often  open  to 
question,  since  the  location  of  ice  observed  from  the  air  is  directly 
dependent  on  the  accuracy  with  which  the  position  of  the  aircraft  can  be 
fixed;  and  concentrations  depend  on  the  accuracy  of  the  human  eye.  The 
Leningrad  Division  of  the  State  Oceanographic  Institute  has  tested  an 
instrument  designed  to  measure  ice  coverage.    In  the  tests,  during  which 
aircraft  position  could  be  accurately  determined,  it  was  noted  that  in- 
dependent visual  observers  not  only  differed  with  observations  made  with 
the  instrument,  but  among  themselves  as  well.   Errors  in  positioning  the 
ice  boundary  were  up  to  15  miles;  errors  in  area  covered  were  up  to  40$ 
[Betin;  et  al,  11].   The  implication  is  that  even  direct  observations 
may  reflect  the  actual  condition  inaccurately. 

Current  ice  reporting  messages  are  not  readily  adaptable  to  digital 
methods.   A  report  must  be  plotted  on  a  chart  before  data  can  be  extracted 
for  program  inputs.   This  is  a  time-consuming  process  and  results  in  a 
great  loss  in  the  speed  advantage  realized  by  a  machine  method. 

Programming  language. 

The  program  as  written  is  in  the  programming  language  called  FORTRAN  '60 
This  language  is  simple  and  easy  to  use  and  has  speeded  the  production  of 
the  program.   However,  it  has  some  disadvantage  in  its  use;  for  example, 
all  data  must  be  input  via  punched  card  to  magnetic  tape,  and  computed 
values  cannot  readily  be  retained  for  use  in  new  forecasts.   Also,  since 
a  wind- forecasting  program  is  not  available  in  the  language,  all  wind 
forecasts  must  be  extracted  from  prognostic  charts  prior  to  being  input. 
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All  this  results  in  loss  of  time  in  entering  data  into  the  program;  the 
entire  process  requires  more  than  15  hours  for  a  600  by  750  nautical  -mile 
grid,  provided  the  required  ice  and  weather  charts  are  already  plotted. 

Program  step  size. 

Program  step  size  is  only  a  limitation  in  that  it  affects  other  parts 
of  the  calculation.   We  have  mentioned  briefly  transient  states  of  the  ice 
field;   it  is  possible  that  the  wind  field  might  be  such  as  to  make  the 
transient  state  important  even  in  the  24-hour  steps  now  taken,  and  these 
stages  surely  would  be  important  if  the  steps  were  shortened  to  the  point 

that  steady  state  did  not  apply.   The  step  size  is  important  in  the  com- 

dc 
putation  of  -r-.   As  now  programmed  this  quantity  is  obtained  by  a  for- 
ward differencing  method  which  can  introduce  errors  over  a  large  number 
of  time  steps.   If  the  number  of  time  steps  were  increased,  a  significant 

error  probably  would  result,  necessitating  a  revision  of  the  method  of 

dc 
obtaining  t-  .   A  central-point  differencing  scheme  is  suggested  as  a 

means  of  accomplishing  this  if  program  steps  are  increased  in  number. 
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7.   Results. 

This  section  discusses  only  ICEGRID  as  ICEF0R1  simply  solves  Shuleikin's 
equations.   Due  to  the  time  involved  in  extracting  data  for  the  program  input 
only  three  runs  were  made  for  verification.   This  small  number  of  runs  does 
not  permit  statistical  analysis  of  the  forecast  versus  the  observed  ice 
field.   Analysis  is  further  complicated  by  the  fact  that  observations  often 
cover  only  a  small  portion  of  the  area  of  forecast;  hence  forecasts  were 
made  for  the  few  periods  when  verification  data  were  available.   Despite 
these  disadvantages,  several  items  of  significance  are  shown  by  the  results. 
The  effects  of  permanent  currents,  for  example,  are  clearly  demonstrated. 
First  to  be  considered,  however,  is  the  degree  to  which  we  can  depend  on 
the  program  as  a  forecast  of  the  wind-drift  component  of  the  sea  ice  motion. 

In  the  analysis  of  results  it  must  be  remembered  that  only  the  wind- 
induced  component  of  drift  is  predicted   The  immediate  conclusion  under 
this  condition  is  that  ordinarily  one  would  not  expect  the  forecast  to 
conform  exactly  with  the  observed  ice  conditions.   Such  a  conformance 
would  result  only  if  the  drift  due  to  permanent  currents  and  heat  budget 
effects  were  negligible. 

Using  coefficients  of  friction  obtained  from  ICEF0R1,  ICEGRID  produces 

almost  identically  the  drift  which  one  obtains  by  hand  using  the  curves 

of  the  U.  S.  N.  Oceanographic  Office  [4]  for  Baffin  Bay  (these  were  for 

k  =  .003  and  k   =  .013).   Thus,  one  is  assured  that  for  areas  of  10/ 10 
a  w 

concentration  forecasts  of  wind  drift  will  be  at  least  as  reliable  as 
those  obtained  by  hand.   Considering  that  we  alter  Shuleikin's  equations 
for  concentration,  we  may  extend  this  reasoning  to  say  that  forecasts 
of  wind  drift  ought  to  be  better  than  those  obtained  by  hand  since  hand 
forecasts  are  made  using  curves  for  10/10  concentration,  without  modification. 
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While  no  hand  forecast  has  been  made,  for  such  comparison  one  machine  fore- 
cast was  done  considering  all  concentrations  to  be  10/10  for  drift  calcu- 
lations; a  significant  drift  difference  was  noted  from  a  forecast  modified 
for  concentration,  using  the  same  input  data.   The  forecast  considering 
all  concentrations  was  better  for  all  points  in  the  grid  than  the  one 
assuming  ice  concentration  to  be  10/ 10  everywhere , 

To  further  establish  the  worth  of  ICEGRID  in  predicting  wind  drift, 
an  examination  of  the  machine  forecast  versus  a  forecast  based  on  persis- 
tence of  the  input  field  was  made.   Figures  3,  4  and  5  show  the  forecast 
boundary  versus  input  and  observed  boundaries  for  each  of  the  three  fore- 
casts made.   In  each  forecast  the  machine  result  was  better  than  persistence, 
in  that  it  gave  results  which  were  generally  more  like  the  observed  ice 
field  so  far  as  ice  boundaries  were  concerned.   The  movement  was  generally 
in  the  right  direction.   Comparing  forecast  concentrations  to  persistence 
of  concentration  also  shows  the  forecast  to  be  better  since  it  gives,  at 
least,  a  hint  as  to  the  direction  of  change  in  concentration.   Figures  6 
and  7  contrast  an  input  concentration  to  its  five-day  forecast. 

In  comparing  the  machine  forecast  with  the  observations  at  verifying 
time  one  finds,  in  general,  a  considerable  discrepancy  in  the  location  of 
ice  boundaries,  say  five  to  30  miles  per  five  days,  and  sometimes  in  con- 
centrations, as  much  as  six  tenths,  which  must  arise  out  of  effects  other 
than  pure  wind  drift.   In  most  instances,  the  forecast  drift  is  of  the  right 
sense,  but  it  is  generally  insufficient  to  produce  the  observed  results. 
A  full  analysis  of  the  forecast  results  versus  observation  requires,  then, 
a  consideration  of  factors  ignored  in  the  development  of  the  method   Most 
of  these  have  been  previously  discussed  (see  section  6)  and  will,  therefore, 
be  mentioned  here  only  as  they  seem  to  affect  the  validity  of  a  forecast 

made  using  only  wind  drift. 
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FIGURE  5:  Limits  of  ice-inpu^forecast&observed-for  June 

11  &16, 1953, showing  use  of  a  geographical  overlay. 
Forecast  area  =600x750  nautical  miles. 
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FIGURE  7 .  FIVE  DAY  Forecast  (concentration)  May  17,1959 
Area  of  forecast  =  600x750   nautical  miles. 
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The  effect  of  permanent  currents  is  the  most  apparent  of  all  effects 
exterior  to  the  wind.   Forecasts  generally  show  the  observed  ice  boundary 
to  vary  from  that  forecast  by  an  increment  which,  for  Baffin  Bay,  appears 
to  conform  well  with  the  permanent  current  believed  to  exist.   It  is  con- 
ceivable that  for  extremely  light  winds  the  effects  of  the  permenant 
current  would  dominate  leading  to  a  drift  entirely  different  from  that 
forecast .   In  running  forecasts,  it  was  noticed  that  in  central  Baffin 
Bay  around  67  N  latitude  forecasts  conformed  rather  well  with  observations, 
indicating  a  lack  of  permanent  current  there,  or  at  least  a  current  parallel 
to  the  ice  boundary.   The  study  of  permanent  currents  through  use  of  this 
program  is  suggested  as  a  possible  research  problem. 

A  second  noticable  effect  is  that  due  to  melting,  since  the  forecasts 
were  all  made  during  the  melting  season.   This  is  seen  most  easily  through 
comparing  observed  concentrations  with  those  forecast.   Large  decreases 
beyond  those  forecast  are  sometimes  observed   Such  changes   are  likely 
due  to  melting  and  not  to  divergence  or  advection  alone,  as  no  reasonable 
amount  of  drift  could  account  for  them. 

In  a  previous  section  data  errors  were  discussed.   These  contribute 
to  apparent  errors  in  the  forecast.   Since  there  is,  at  present,  no  method 
available  to  exclude  such  data  errors,  one  must  use  all  data  assuming  it 
to  be  correct.   Fortunately,  with  regard  to  ice  concentration,  it  is  the 
navigability  of  an  area  which  is  important  and  not  the  absolute  ice  con- 
centration.  Thus,  some  error  in  forecast  concentration  may  be  tolorated. 
The  forecast  concentration  has  value  if  it  produces  an  indication  as  to 
the  navigability  of  the  frozen  sea.   Analysis  of  program  outputs  indicate 
the  method  to  be  capable  of  such  indication.   In  the  forecasts  produced 
it  has  always  been  possible  to  determine  rather  accurately  areas  of  de- 
creasing concentration*  hence,  one  is  able  to  determine  the  areas  most  likely 
to  be  navigable.  9Q 


Analysis   of   forecasts    indicates    the    input  wind    fields    to  be   poorly 
defined  with   respect   to  divergence   and  curl   of   the  wind   on  the   scales 
desired.      Thus,    small-scale    features    in   the    ice    field   are   not   shown   or  are 
smoothed   out   over   the    period   of   forecast.      These   small    features   are   often 
transient  and  highly   important    in   the    local  navigation  problem       This   raises 
a   question  as    to  the   effect   of  small-scale    features    in  the  wind    field 
on  the  movement   of  sea    ice.      Browne   and   Crary    [9]    indicate   that   the    ice 
fields  may  move  with   the  average  wind   over  a   60-nautical-mile   square 
surrounding   the    location,    but  are   not  able    to  determine   conclusively 
this   relation,    as  weather  maps    fail   to  supply  adequate  wind   information 
on  this   scale;  and,   moreover,    all      sea    ice    is    not    in  the    form  of   ice 
fields.      It  may  be,    then,    that  wind   forecasts    to   include   smaller- scale 
features    in   the  wind   field  will  be   necessary  before    the   smaller- scale 
features    of   ice  drift  will  become  apparent.      Since  weather   reports    in 
the  Arctic  areas   are    few   in  number  such  a   definition  of  the  wind   field 
may  never  be   possible. 

The   results    obtained  make    it  apparent   that   areas    of  Baffin  Bay  exist 
wherein  nearly  all   the   drift   of   ice   is    that   due   to  wind.      This   seems    par- 
ticularly true   in  the   central  area  near  and   just   north   of   the   Davis   Strait. 
In  such  an  area   the  wind-drift   forecast   should   constitute   a   total   forecast 
if  no  heat   budget   effects   are   present. 

The   few  forecasts   made   using   program  ICEGRID   indicate   that   the  method 
produces   good   results    for  wind  drift  and   that   it    is    possible    to  expand    it 
to   include   other  effects. 
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8.   Recommendations. 

A  review  of  the  preced  ing  sections  invites  consideration  of  ways  of 
improving  the  program.   One  is  struck  by  the  need  for  better  methods  of 
obtaining  and  reporting  data.   Additionally,  results  obtained  point  out 
the  necessity  of  expanding  the  method  to  provide  a  complete  ice  forecast. 

A  study  of  methods  of  observation  and  reporting  of  data  has  not  been 
attempted  here  and  might  be  the  subject  of  separate  investigation.   A 
few  possible  aids  to  improvement  of  data  quality  and  availability  are 
worthy  of  note,  however.   Stereographic  aerial  photography  and  TIROS 
(Weather  Satallite)  photographs  are  suggested  as  providing  more  complete 
observational  data,   Stereo  photographs  would  give  detail  as  to  the  dimensions 
of  features  in  the  ice  such  as  hummocks,  pressure  ridges,  and  even  an  estimate 
of  thickness  where  ice  and  open  water  lie  adjacent.   Large  scale  concentra- 
tion data  should  be  available  from  TIROS  photographs  where  cloud  cover  does 
not  obscure  the  surface.   Scanning  of  photographs  might  be  done  electronically, 
obtaining  directly  concentration  inputs  for  each  grid  point. 

Even  with  better  methods  of  obtaining  data,  time  required  to  input 

that  data  to  ICEGRID  will  remain  large  if  the  reporting  format  is  not 

changed.   The  present  format  presents  location,  amounts  and  characteristics 

of  the  ice  verbally,  requiring  the  data  to  be  plotted,  extracted  for  each 

grid  point  and  punched  onto  data  cards  prior  to  computer  input.   A  suitable 

format,  from  which  a  direct  input  to  the  computer  could  be  obtained,  might 

be  similar  to  that  of  the  World  Meteorological  Organization  (WMO)  used  to 

report  weather.   Such  a  report  format  is  suggested  in  figure  8  and  could 

be  included  in  the  WMO  code   This  format  would  report  total  concentration, 

concentrations  by  size  and  boundary  location  points  for  the  concentration, 

all  reported  digitally.   Additional  plain  language  data  could  be  included 

at  the  end  of  such  a  message.   The  code  groups  could  be  repeated  for  each 
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C 'C  hi  ■  rij  n,    LaLaL0L0t    LaLaL0L0t 

C  C  .'  concentration  in   tenths 
Pl1       :  tenths  of  slush,  brash,  block 
Pl2      .'  tenths  of   small  to  medium  floes 
n3     ;  tenths  of  large   f loes  &  field 
LaLaMatitude  in   whole   degrees 

L0L0*  longitude  in  whole  degrees  as 
in  World  Meteorology  code 

t        !  tenths   of   degrees  indicator  as 
below: 
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9 

FIGURE  8:    Suggested  ice    reporting  format 
to  be  used   with  the    World 
Meteorology    code. 
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"concentration  area" of  ice  observed.  The  boundary  between  ice  and  open 
water  could  be  described  by  similar  latitude-longitude  groups  following 
a  code  group,  LIMIT. 

We  must  also  seek  to  provide  a  means  for  machine  retention  of  com- 
puted data  and  a  wind -forecasting  program.   Fortunately,  the  U.  S.  Navy 
Fleet  Numerical  Weather  Facility  has  available  a  program  computing  a  wind 
prognosis.   This  suggests  rewriting  ICEGRID  in  SCRAP,  the  programming 
language   used  by  FNWF,  obtaining  immediately  a  machine  forecast  of  winds 
and,  at  the  same  time,  the  ability  to  retain  computed  data  as  continuity 
for  new  forecasts. 

Extension  of  the  program  to  include  the  effects  of  permanent  currents, 
heat  budget  effects,  land  and  fast=ice  effects  on  the  drift  vector   (in 
short,  a  complete  forecast)is  also  suggested.   Such  expansion  would  require 
considerable  time  and  effort  on  programming  and  on  obtaining  the  necessary 
topography  and  oceanographic  data,   The  production  of  an  efficient  and 
useful  forecast  method  dictates  the  need  for  this  effort,  however,  and  en- 
courages that  time  be  spent  on  such  expansions. 
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9.   Conclusions. 

In  the  foregoing  sections  the  programs  presented  here  have  been 
examined  and  some  of  the  factors  foreign  to  pure  wind  drift  of  sea  ice, 
which  one  ought  to  consider  in  a  complete  forecast,  have  been  discussed. 

The  programs  developed  here  are  concluded  to  be  superior  to  hand  fore- 
casts on  speed  and  on  accuracy,  due  to  the  modification  made  for  concentra- 
tion and  the  ability  to  retain  computed  values  to  provide  continuity 
to  the  forecasts.  With  the  inclusion  of  the  effects  of  permanent  currents 
and  heat  budget,  one  may  expect  to  obtain  a  complete  and  useful  ice-drift 
forecast,  provided  good  inputs  are  available. 

Further,  it  is  concluded  that  the  ability  to  obtaia  values  rapidly  for 
the  coefficients  of  friction,  k  and  k  ,  leads  to  the  programs1  great  value 
in  research  work. 
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APPENDIX  A 
PROGRAM  ICEFOR1 
Description. 

ICEFOR1  utilizes  Shuleikin's  equations,  without  modification,  to 
compute  the  wind  drift  of  sea  ice.   Provision  is  made,  however,  for  entry 
of  coefficients  of  friction  suitable  to  the  area  of  computation.   Drift 
vectors  are  computed  for  one  or  more  locations  provided  the  primary  in- 
puts of  latitude,  ice  thickness,  wind  speed  and  wind  direction  are 
entered.   Secondary  inputs  of  longitude  and  date/time  group  may  also  be 
entered  for  reference  purposes,  but  these  are  not  required  in  any  cal- 
culation. 

This  program  has  its  primary  use  as  a  research  tool„   Given  a  set  of 
observations  of  ice  drift,  corrected  for  permanent  currents,  and  of  corres- 
ponding wind  velocities,  one  may  find  the  Shuleikin  curve  with  the  pair 

of  coefficients  of  friction,  k  and  k  ,  which  best  fit  the  data.   Having 

'aw  ° 

obtained  the  constants,  one  may  obtain  ice-drift  vectors  to  be  used  in 
forecasting . 

Table  A-I  defines  the  variables  used  in  the  program  language.   Where 
variables  are  equivalent  to  those  in  Shuleikin's  equations  symbols  are 
included.   The  complete  program  is  listed  following  this  discussion. 
Data  input. 

Prior  to  use  of  the  program,  it  is  necessary  to  assign  values  to  the 
constants  AKA  and  WKW  in  the  subprogram,  DRIFOR.   The  values  assigned  these 
constants  in  Figure  A-l  are  for  Baffin  Bay.   The  values  assigned  to  RHOAIR 
and  DSUBW  may  also  be  changed,  but  this  is  not  usually  required. 

Data  input  is  accomplished  using  IBM  cards  with  punched  values  of 

variables  entered  on  them„   The  first  card  specifies  latitude,  longitude, 

date/time  group,  ice  thickness  and  wind  velocity,  assumed  constant  for  three 
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hours.   The  second  and  subsequent  cards  are  available  for  specification 
of  additional  three-hourly  wind  data  at  the  previous  data  point.   Follow- 
ing the  last  card  specifying  wind  data  is  a  card  on  which  wind  speed  is 
entered  as  5000.   This  signals  the  computer  to  proceed  to  the  next  card 
specifying  the  complete  data,  i.e. :  another  "first  card".   When  all  data 
have  been  read,  a  card  which  specified  latitude  as  300  signals  the  computer 
to  read  no  more  data  cards „ 
Output  of  results 

The  program  prints  output  in  columns,  titled  as  follows: 

Result  Units 

latitude  degrees 

longitude  degrees 

date-time  day-hour 

ice  speed  knots 

ice  direction  degrees 

distance/3  hrs „  nautical  miles/day 

psi  radians 
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TABLE  A-I 

Definitions  of  variables  in  ICEF0R1 

PLAT  -   latitude 

PLONG  =  longitude 

DTG  =  data  and  time  group 

THICK  =  ice  thickness 

VELWIND  =  wind  speed  in  knots 

DIR  WIND  =  wind  direction 

VELICE  =  ice-drift  speed  in  knots 

DIRICE  =  ice-drift  direction 

DIST3  =  3- hour  drift  distance  in  nautical  mile 

AKA  =  coefficient  of  friction  between  air  and  ice,  k 

WKW  =  coefficient  of  friction  between  ice  and  water,  k 

RHOAIR  =  air  density,  6a 

DSUBW  =  water  density,   6 

ENN  =  N  W 

OMEGAB  =  u)  sin  to   ■  a> 

XI  =  § 

ALFA  =  a 

BETA  -  p 

SQMASS  =  M 

BEE  »  B 

PSIO  =  an  angle,  used  in  computation 

PSI  =  angle  between  wind  and  ice-drift  direction,  i|r 

WICE  =  ice-drift  speed  in  meters  per  second,  u 
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APPENDIX  B 
PROGRAM  ICEGRID 
Description. 

ICEGRID  computes  the  wind  drift  of  sea  ice  using  equations  after 
Shuleikin,  as  modified  for  concentration  less  than  10/10,  and  concentra- 
tion changes  according  to  the  equations  of  Simpson.   Drift  vectors  are 
computed  over  a  546=point  array  using  latitude,  ice  thickness,  ice  con- 
centration, wind  speed  and  wind  direction  as  inputs.   The  purpose  of  the 
program  is  to  obtain  a  five-day  forecast  of  the  wind  drift  of  sea  ice. 

Table  B-I  defines  the  variables  used  in  the  program  language.   Where 
variables  are  equivalent  to  those  of  Shuleikin  or  Simpson,  symbols  are 
included „   The  complete  program  is  listed  following  this  discussion. 
Data  input 

Data  input  is  accomplished  using  IBM  cards  with  the.  values  of  the 
variables  punched  on  them0   Four  constants  must  also  be  entered  into  the 
program;  these  are  also  input  via  IBM  card. 

The  first  data  card  specifies  the  constants,  SCALE,  which  is  the 
scale  of  the  grid,  and  BLAT,  the  base  latitude.   A  scale  of  one  provides 
grid  squares  of  60  nautical  miles  on  a  side;  a  scale  of  two  gives  grid 
squares  30  nautical  miles  on  a  side;  and  so  on.   BLAT  establishes  the 
northern-most  boundary  of  the  grid  and  is  always  the  latitude  of  that 
boundary  plus  the  inverse  of  the  scale  factor.   The  second  data  card 
specifies  the  values  assigned  to  the  coefficients  of  friction,  AKA  and  WKW. 
The  fourteen  data  cards  following  the  two  just  defined  specify  the  concen- 
tration of  ice  at  all  of  the  546  grid  points,  values  determined  by 
averaging  observations  over  the  four  adjacent  grid  squares.   Land  points 
are  entered  as  2C0.   Following  these  cards  are  fourteen  cards  specifying 

ice  thickness  over  the  grid,  taken  to  be  the  mean  climatological  thickness 
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for  the  period  of  forecast,,   Thickness  is  entered  in  feet  and  is  obtained 
in  the  same  manner  as  concentration  for  each  grid  point.   Fourteen  more 
cards  specify  wind  speed  in  knots,  obtained  from  prognostic  charts,  for 
each  grid  point  by  averaging  over  the  four  adjacent  grid  squares.   Fol- 
lowing this  group  of  cards  is  a  group  of  21  cards  specifying  the  direction 
toward  which  the  wind  is  blowing. 

The  remainder  of  the  input  cards  consists  of  groups  of  14  and  21  cards 
respectively  for  wind  speed  and  direction  for  each  of  the  remaining  fore- 
cast days. 
Output  of  results 

Results  are  printed  in  grid  form  with  the  outer  boundary  points  of 
the  grid  omitted.   The  input  concentration  grid  is  reproduced  as  the  first 
output.   This  is  followed  by  grids  of  concentration,  direction  of  movement 
and  distance  for  each  of  the  five  forecast  days. 
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TABLE  B-I 


Definitions  of  variables  in  ICEGRID 


CONCICE  -  ice  concentration,  c 

VELWIND  =  wind  speed  in  knots .  V* 

DIRWIND  ■  wind  direction  +180  ,  in  degrees 

THICK   =  ice  thickness,  enters  into  computation  of  mass 

TEMPICE  =  a  temporary  storage  for  new  concentration 

VELX    =  x-component  of  ice  speedy  u,  in  knots 

VELX    =  y-component  of  ice  speed,  r,  in  knots 

DIST  24  =  distance  of  ice  drift  in  24  hours,  in  nautical  miles 


=  scale  factor,  S, determines  grid  size 

=  base  latitude  permits  proper  orientation  of  grid  in  latitude 

=  latitude,  degrees,  0 

=  ice  drift  speed  in  knots 

=  ice  drift  direction,  in  degrees 

=  change  in  concentration  with  time,  dc/dt 


=  coefficient  of  friction  between  air  and  water,  k, 
=  coefficient  of  friction  between  ice  and  water,  k 


AKA 

WKW  =  coefficient  of  friction  between  ice  and  water,  k~ 

RHOAIR  =  air  density,  6a 

DSUBW  =  water  density,  6 

OMEGAB  =  go  sin  0,  SB 

XI  =  § 

ENN  =  N 

ALFA  =  a 

BETA  -  0 

SQMASS  ■  m 

BEE  -  B 

PSIO  =  an  angle  ,  used  for  computation 

PSI  =  angle  between  wind  and  ice  direction,  i|r 

CONFAC  =  factor  modifying  for  concentration 

VVICE  =  ice  drift  speed  in  meters  per  second,  u 

C  =  constant  used  to  modify  outputs  of  i|r  and  u,  if  consistently 
low  or  high,  normally  set  to  one 
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